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Abstract. In this paper we investigate a discrete approximation in time and 
in space of a Hilbert space valued stochastic process {u(t)}tg[o,T] satisfying a 
stochastic linear evolution equation with a positive-type memory term driven 
by an additive Gaussian noise. The equation can be written in an abstract 
form as 

du +(^j b(t - s)Au(s)ds S J dt = dWQ , t £ [0,T]; u(0)=u eH, 

where is a Q- Wiener process on H = L (T>) and where the main example 
of b we consider is given by 

6(t) =t^- 1 /r(/9), 0</3<l. 

We let A be an unbounded linear self-adjoint positive operator on H and we 
further assume that there exist a > such that A~ a has finite trace and that 
Q is bounded from H into D(A K ) for some real K with a — -gq-j- < n < a. 

The discretization is achieved via an implicit Euler scheme and a Laplace 
transform convolution quadrature in time (parameter At = T/n), and a stan- 
dard continuous finite element method in space (parameter h). Let u n ^ be 
the discrete solution at T = nAt. We show that 

(E\\u n , h -n(T)|| 2 ) 1/2 =0{h v + M">), 

for any j < (1 — (J3 + l)(c? — k))/2 and v < -o^f — a + k. 



1. Introduction 

Let T> be a bounded domain in R d , d € N, and let a be a real- valued stochastic 
process solution of the equation formally written as 

(1-1) j\(t-s)Au(x ) s)ds = i(x > t), (x,t)eVx(0,T], 

together with the initial condition u(x,0) = Uq(x), x G T>, and boundary condition 
u\dx> = 0. Here, £ is a zero mean real valued Gaussian noise and the time kernel 
b is assumed to be real- valued and of positive type; i.e., that for any T > 0, the 
kernel b belongs to i 1 (0, T) and for any continuous function / on [0, T] the following 
inequality holds: 

(1.2) f T f b(t-s)f(s)f(t)dsdt>0. 

Jo Jo 
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The deterministic version of such problems can be used to model viscoelasticity 
or heat conduction in materials with memory (see |13j for references). When b is 
smooth, these equations exhibit a hyperbolic behaviour, whereas if b has a weak 
singularity at t — (for example a Riesz potential), they exhibit certain parabolic 
features. In particular, when 

(1.3) b(t) = tP- 1 /T(p), 0<p<l, 

the homogeneous deterministic equation has a smoothing property which corre- 
spond to the inequality 

(1.4) ||« (m) (*)llff*(») < Ct-^ r - m \\u \\ L2{v) , 

where |r| < 1 if m > 1 and where 0<r<lifm = 0, but with no further smoothing 
in the spacial variables (sec e.g. [131 Theorem 5.5]). The framework of this paper 
allows for slightly more general kernels but with similar smoothing effects and, in 
particular, they are of positive type. Hence, together with the positivity of the 
operator —A, the deterministic equation will remain parabolic in character. 

Next we introduce an abstract framework to describe the noise and equation (11.11) 
more precisely. Let Q be a self-adjoint, nonnegative linear operator on H := L 2 (D) 
and be a Wiener process in H with covariance operator Q (or, simply, Q- 
Wiener process). We set A = -A, D(A) = H 2 {V) n H&(V) and H = L 2 {V). Then 
A can be seen as an unbounded linear operator on H with domain D(A). For 
b given in ()1.3[) and under reasonable assumptions on 3D, our main assumption 
concerning Q in (11.11) is that A K Q defines a bounded operator on L 2 (D) with 
d/2- 1/08 + 1) < k < d/2. 

If we write u(t) = u(-,t), considered as a ii- valued stochastic process, then (jl.ljl 
can be rewritten in the abstract Ito form as 

(1.5) du(t) + (J b(t- s)Au(s)ds^j dt = dW Q (t), t£(0,T], 

with initial condition u(0) = uq G H. 

While the literature on numerical methods for deterministic infinite dimensional 
Volterra equations is abundant (sec, for example, [TJ [6J [12j [T3l HO], which is a 
very incomplete list), the numerical analysis of stochastic Volterra equations is 
completely missing. We are only aware of [7] where an algorithm is described 
and numerical experiments are performed with no error analysis given. We will 
consider a numerical approximation of (|1.5[) by means of an Euler scheme and a 
Laplace transform convolution quadrature in time together with a finite element 
method in space. Let n > 1 an integer, At = T/n and tk = k At, k = 0, . . . , n. Let 
also {Vh}h>o t> c a family of finite dimensional subspaces of D(A X / 2 ) = Hq(V). For 
each 1 < k < n, we seek for an approximation of u(tk) in Vh by Uk,h defined by the 
following induction: 

k 

(1.6) (uk,h - u k -i,h,Vh) + At22wk-j[Auj,h,Vh) = ^ r ^~i{Q 1/2 Xk,Vh), k > 1, 

for any vh S Vh, where \J At \k is the noise increment and where (-,-) is the inner 
product of H . The approximation of the convolution term in (|1.5j) is achieved via 
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a quadrature rule such that for any continuous function / on [0, T], 



«»k-jf(tj)~ / b(t k -s)f(s)ds = (b-kf)(t k ). 



Then, the approximation of 6*/ on the time grid t k , k = 0, . . . , n, is obtained from 
a discrete convolution with the values of / on the same grid. Before going into 
details, let us point out that not any quadrature rule can be chosen. In particular, 
it will be important for the chosen quadrature to satisfy a discrete analogue of (|1.2I) . 

In order to understand the specific quadrature rule used in this paper, we will 
take the example of the Riesz kernel (|1.3[) . Let us note that in this case the Laplace 
transform of b is z~@ and the term b * Ait in (|1.1|) can be seen as the fractional 
integral (d/dt)~"(Av). Then, the idea is to use the same Euler approximation 
of d/dt in both terms on the left hand side of (|1.1[) . Since the discrete Laplace 
transform of the implicit Euler scheme is (1 — z)/At, one chooses the quadrature 
weights to have discrete Laplace transform ((1 — z)/At)~@. 

Such a convolution quadrature has been introduced in [9j HOj . It was motivated 
by the fact that the main properties of the solution of the homogeneous problem, 
like stability, existence, or regularity, are largely determined by the distribution of 
the frequencies of the kernel (by means of its Fourier or Laplace transform) , espe- 
cially when the kernel has weak singularities or when it exhibits different behaviour 
at different time scales. Since, by construction, the discrete Laplace transform of 
the quadrature kernel is closely related to the Laplace transform of the continuous 
kernel, it is thus possible to carry over frequency domain conditions from the con- 
tinuous problem to the discretization and thereby obtain stable approximations. 
Moreover, this kind of quadrature rule inherits the rate of approximation from the 
time integrator of d/dt. In the context of stochastic PDEs, we think that it is 
important to make sure that the deterministic part of the scheme is stable and that 
the perturbations are due to the noise only. 

Although the analysis in the present paper allows for kernels slightly more general 
than (|1.3p . we follow the same idea: the convolution quadrature weights {co k } in 
(|1.6j) will be defined by means of the Laplace transform of the kernel b. Therefore, 
we choose the quadrature coefficients to have generating function b((l — z)/At) 
where b denotes the Laplace transform of b; that is, 



We will not focus here on practical algorithms for the computations of the quadra- 
ture weights and we refer the reader to, for example, |10j . 

While precise conditions on the kernel b are postponed to Sections [2] and HI we 
can already state our main result, Theorem l5.ll with the above notations in the case 
of the specific kernel (| 1 . 3[) when I? is a convex polygonal domain using continuous, 
piecewise linear finite elements. We shall prove a (strong) error estimate of the 
form 



where 7 < (1 - (/3 + l)(d/2 - k))/2 and u < 1/(0 + 1)- d/2 + k. Let us note that 
we recover the known order of convergence for the heat equation (see [5J [HI [3T]) 
when j3 — > 0. 



(1.7) 




(E(||u ni/l 



u(T)\\ 2 ) 1/2 <C(Af +h v ), 
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The paper is organized as follows. In Section [2] we introduce notations, recall 
some basic preliminary results, and state our main assumptions on A, Q and b. 
We note that Assumptions (|2.8p - (|2.9[) on A and Q could be replaced by a single, 
somewhat sharper, assumption as discussed in Remarks 12.81 13.51 14.71 and 15.21 It 
is, however, harder to check in most cases. In Section [3] we study the space semi- 
discretization of and strong error estimates are derived for smooth initial 
data under minimal regularity assumptions (Assumption [l| on b. In Section |4] we 
prove strong error estimates for the time semi-discrete scheme with non-smooth 
initial data. One of the key results in this direction is Theorem 14. 1[ where we 
prove a general P-stability result on Lubich's convolution quadrature based on the 
Backward Euler method for deterministic Volterra equations. Interestingly, this 
stability result implies (Corollary I4.2|) that the time-discrete scheme exhibits the 
same smoothing effect in time as the solution under Assumption [T] on b. However, 
in order to obtain optimal convergence rates for the stochastic problem we need to 
put a further regularity restriction on b in Subsection 14.21 Assumption [2] which is 
in fact common in the deterministic literature for nonsmooth initial data. Indeed, 
Assumption [2] implies that the deterministic equation has an analytic resolvent 
family while Assumption[T]only implies that the deterministic equation is parabolic. 
Unlike for equations with no memory term, these two notions are not equivalent (See 
|17l Chapter 1 , Section 3] ) . As far as we know the derivation of nonsmooth initial 
data estimates using only parabolicity (Assumption [1} remains an open problem. 
Finally, in the last section, we gather the results from the preceding sections and 
consider the fully discrete scheme. 

2. Notations and preliminairies 

Let (X, || • || jsc ) and (Y, || • ||y) be two Banach spaces and let B(X,Y) denote 
the space of bounded linear operators from X into Y endowed with the norm 
ll-5||g(A-,Y) = sup^g^ ||Ba;||y/||a;||x- When X = Y, we use the shorter notation 
B(X) for B(X, X). If A is a Banach space and I is an interval in K then, L P (I, X), 
1 < p < oo, denotes the space of functions / : I — > X which are measurable and 
t — > \\f(t)\\ p is integrable on /, equipped with the usual norm. If p = oo then 
i°°(7,A), denotes the space of functions / : I X which are measurable and 
(4 ||/(i)|| is essentially bounded on I endowed with the usual norm. 

Throughout this paper, H denotes a real separable Hilbert space with inner 
product (•, •) and associated norm || ■ || . We consider the stochastic Volterra equation 
given in the abstract Ito form as 

(2.1) du + ^ b(t - s)Au(s) ds^j dt = dW Q , t G (0,T]; u(0) = u G H, 

where the process {w(£)}te[o.T] is a -ff-valued stochastic process, A is a densely 
defined, nonnegative self-adjoint unbounded operator on H with compact inverse, 
and W Q is a Q- Wiener process in H on a given probability space (f2, T, P). The 
weak solution of (|2.ip is a mean-square continuous if- valued process satisfying 

(u(t),r))+ [ f b(r- s)(u(s),A*r)) dsdr = (u ,r)) + f (r), dW Q (s)) , 
Jo Jo Jo 

for all j] G D(A*) almost surely for all t e [0,T]. 

It is well known that such assumptions on A implies the existence of a sequence 

of nondecreasing positive real numbers {\k}k>i and an orthonormal basis {ek}k>i 
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of H such that 

(2.2) Ae k = A fe e fe , lim X k = +00. 

k— > + oo 

We define classically, by means of the spectral decomposition of A, the domains 
D(A S ) of fractional powers s £ R of A and we set 

\\v\\ s = \\A s ^v\\, veD(A s / 2 ). 

Remark 2.1. We note that since A is nonnegative self-adjoint, —A generates an 
analytic contraction semigroup on H. Moreover, for any 9 < n, there exists Mg > 1 
such that the following holds: 

M 

\\(zl + Ay l \\ B(H) < — , for any z e So, 

where L 9 = {z£ C\{0}, |arg(z)| < 9}. 

Let Ci(H) denote the set of nuclear operators from H to H; that is, T <G C\(H) 
if there are sequences {oj}, {&,} C if with Y^jLi ll a jllll^jll < 00 an d sucn that 

OO 

T.t = \^(x, bj)aj, x e if. 

3 = 1 

Sometimes these operators are referred to as trace class operators. For T G £i(if) 
we define Tr(T), the trace of T, by 

Tr(T) = J2(Be n ,e n ), 

n=l 

where {e„} is an orthonormal basis of if. This definition turns out to be inde- 
pendent of the choice basis. Furthermore, if L G £i(if) and M G £>(if), then 
LM, ML G £i(if) and 

(2.3) Tr(LM) = Tt(ML). 
If L is also symmetric and nonnegative, then 

(2.4) Tr(ZM) < Tr(L)||M|| B(if) . 

Hilbcrt-Schmidt operators play also an important role in this paper. An operator 
L G B(if) is Hilbert-Schmidt if L*L G £i(ff) or, equivalently, ii* G £i(if). We 
denote by £2 (if) the space of such operators. It is a Hilbert space under the norm 

(2.5) \\L\\c 2{ H) = (Tr(LU)) 1 / 2 = (Tr(LL*)) 1/2 . 

Our analysis will also use the Laplace transform. Let / : M+ — > if be subexpo- 
nential; i.e., that for any e > the function t H> f(t)e~ £t belongs to L 1 (IR+,iJ). 
We define the Laplace transform of / : C+ — > if by 

r + ao 

f(z) = / f(t)e- zt dt, Rez > 0, 
Jo 

where we have used the same notation if for the complexification of if. We denote 
by * the Laplace convolution product on [0, t] of two locally integrable subexpo- 
nential functions f,g G Lj oc (M. + ,H) defined as 

(/*ff)(*)= f f(t-s)g(s)ds. 
Jo 
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It is well known that / * g £ Lj (M.+ , H) is subexponential and 

fkg (_) = f(z)g(z), Rcz>0. 

2.1. Main assumptions. Next we state the main assumptions concerning the 
kernel b and the operators A and Q, which will be used throughout this paper. 

Regarding b, first note that property (|1.2[) can be characterized by means of the 
Laplace transform b of b. It is equivalent to say that Rc(&(A)) > for any Re A > 
(see [H] or |TTJ page 38]). Now it is clear that the positivity property ()1.2|) is 
not sufficient, in general, to ensure smoothing effects like (| 1 -4[) when working with 
kernels that are more general than (|1.3p . This is why, following [3] and [14] , we will 
impose stronger conditions on b. 

Assumption 1. The kernel ^ b £ Lj oc (W+), is 3-monotone; that is, b, —b are 
nonnegative, nonincreasing, convex, and lim.t_j.oo b(t) = 0. Furthermore, 

(2.6) p := 1 + -sup{|arg6(A)|, Re A > 0} £ (1,2). 

In the special case of the Riesz kernel given in (|1.3[) one can easily show that 
p = 1 + /3. From now on we set /? = p — 1 with p defined by (|2.6[) . 

Remark 2.2. It follows from jTTJ Proposition 3.10] that for 3-monotone and locally 
intcgrable kernels b, condition (|2.6[) is equivalent to 

7 fn sb(s) ds 

2.7 lim 7° . < +00. 

^f*-sb(s)ds 

Also note that, by (|2.6p . we have that Re(6(A)) > for Re A > and hence b 
satisfies (|1.2[) . 

For A and Q we suppose that there exists numbers a > and n € K such that 

(2.8) Tr(A- a ) < +00, 

(2.9) A K Q e B(H), a- - < k < a. 

P 

2.2. The nonhomogeneous deterministic problem. Given / £ L 1 ([0, T]; 7?), 
Assumption [T] together with the fact that A is positive and sclf-adjoint implies that 
the deterministic problem, 

(2.10) u(t)+ f b(t- s)Au(s)ds = f(t), ie(0,T], _(0) = u £ 

is well posed for all T > 0. Indeed, there exists a resolvent family {S(t)} t >o £ S(if) 
which is strongly continuous for t > 0, diffcrcntiablc for < > and uniformly 
bounded by 1, see [17l Corollary 1.2 and Corollary 3.3]. The unique mild solution 
of (|2.10j) is given by the following variation of parameter formula |17l Proposition 

1-2] 

u{t) = S(t)u + I S(t-s)f(s)ds, te[0,T]. 
Jo 

Remark 2.3. The positivity of the kernel b defined in (|1.2[) , together with the 
positivity of the operator A already allows for the construction of a unique solution 
to (|2.10p using an energy argument, see [T7J Corollary 1.2]. Assumption Q] gives 
further integrability and smoothing properties for {S(t)}t>o- 
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Note that such a resolvent family does not satisfy the semi-group property be- 
cause of the non local feature of the memory term in (|2.10j) . Nevertheless, it can 
be written explicitly using the spectral decomposition (|2.2I) of A as 



+00 

(2.11) S(t)v = J2s k (t)(v, 

fc=i 

where the functions Sk{t) are the solutions of the ordinary differential equations 

(2.12) s k (t) + X k [ b(t-s)s k (s)ds = 0, s fc (0) = l. 

Jo 

The next proposition summarizes the main properties of the functions {sk}k>i- 

Proposition 2.4. Suppose that b satisfies Assumption]]] and let p G (1,2) as 
defined in i2.6\ ). Then lim r _>. 00 Sfc(r) = for all k > 1 and there exists Co > such 
that for any k > 1. 

(2.13) IMU~(R+) < 1, 

(2.14) l|4|Ui( R+) < C 0j 

(2.15) P4|Ui( R+ ) < C \- 1/p , 
(2-16) \\s k \\mm + ) < C \; 1/p . 

Proof. Estimate (|2~13)) follows from [HI Corollary 1.2], inequalities (|2~14|) and ([2~T5|) 
can be found in [TH Proposition 6] and estimate (|2.16[) is shown in [31 Lemma 3.1] 
where also the fact liuir^aa Sfc(r) = for all k > 1 is shown in the proof of the 
lemma. □ 

Smoothing effects of the resolvent family {S(t)} t >o when b satisfies Assumption 
[T] can be now easily proved using Proposition [ 



Proposition 2.5. Let b and p as in Proposition \ 2.J\ Then for any t > 0, there 
exist a constants Co, C\ > such that for any < s < 2/ p and < s' < 2, 

(2.17) \\A^S(t)\\ B(H) < C a t~ s P/\ t>0, 

(2.18) \\A-°'/ 2 S(t)\\ B(H) < CmUl,/' 2 - 1 , t>0 - 

Proof. For any 5 G (0,1) and any k > 1, Holder's inequality, (|2.14j) and (|2.15[) 
yields 

+ 00 r-\-oo 

u 5 \s k (u)\du = / u^Skiu^lsKu)] 1 - 6 du 
Jo 

< / u\s'k(u)\ du I / |s' fe (u)|du 







< Co A-^. 

Note, that the previous final estimate also holds for S = 0, 1 by (|2.14[) and (12.151) . 
Then, since Sk{t) = — J t + °° u~ s u s s'k(u) du as lim r ^.oo Sfe(r) = for all fc > 1 by 
Proposition ^. 41 we can conclude that 

(2.19) |* fc (t)| <C Q t- s X~ 5/p , t>0, 5e[0,l]. 
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Thus, for any s G [0, 2/p] and x G H, (f2~T9|) with < 5 = ps/2 < 1 implies 
\\A s ' 2 S{t)x\\ 2 = J2 X kSk(t) 2 (x,ekf <C t-<> s / 2 \\xf, 



fc>i 



which is (|2.17[) . To show (|2.18[) . we use [TTJ Corollary 3.3] which states that under 
Assumption Q] and since belongs to the resolvent set of A, there is M > such 
that 

(2.20) \\S(t)x\\ < Mt^WxW, xeH, t>0. 

On the other hand, we can bound S(t)x for x <G D(A) as follows: 

(2.21) ||^)x|| 2 = ^( Sfe (t)) 2 (^ efc ) 2 

fc>i 

(2.22) =]TA 2 ( f t b(t-a)a h (s)da) (x, e k f < \\bf L1{0 , t) \\Ax\\ 2 , 

k>i \ Jo ' 

where we have used (12.121) and (|2.13j) . Finally, interpolation between (|2.20|) and 
(12~2TT) yields (l2~T8l) . □ 

Remark 2.6. The estimate in (|2.18j) does not provide an optimal rate, in fact it is 
the worst case scenario, as further smoothing may come from ||6||ii(o,t)- The rate 
can be improved if we impose further regularity assumptions on b. Indeed, if in 
addition, b satisfies Assumption [2] from Subsection 14.21 then by (|4.3|) and (|4.8[) it 
follows that 6(A) ~ A 1- '' as A — > oo. Thus, the nonnegativity of b implies that 
II^IIl^o,*) Ci p_1 by a Tauberian theorem for the Laplace transform (see, for 
example, [22j Chapter V, Theorem 4.3]). Therefore, in this case, we get a sharper 
estimate 

\\A- s '/ 2 S(t)\\ B{H) < d t ps '/ 2 -\ t > 0, < s' < 2. 

Nevertheless, the rate in given (|2.18[) is sufficient for our needs when it is used in 
the deterministic error analysis for smooth initial data. 

2.3. The continuous stochastic problem. Next we recall an existence result for 
the problem (|2.1j) and, for the sake of completeness, we indicate a proof (see [H 
Theorem 2.1] and we refer to |18j for more general noise). 

Proposition 2.7. Let A and Q satisfy ()2.8j) - (12.9[) and let b satisfy Assumption^ 
Then there exists an unique H -valued (Gaussian) weak solution u of (|2.1j) given by 
the variation of constants formula 

(2.23) u{t) = S(t)u + [ S(t- s)dW Q (s). 



Furthermore, the stochastic convolution term has a version whose trajectories are 
a.s. 6-Hdlder continuous for any 9 < (1 — p(a — k))/2. 
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Proof. Analogously to [U Theorem 5.4], it is sufficient to show that the stochastic 
convolution is well-defined. By Ito's Isometry, 

2 rt 

\\S(t-s)Q 1/2 \\ 2 c 2 (H } d S 
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S(t- s)dW Q (s) 



f y2\\S(t-s)Q 1 / 2 e i fds 
Jo i>i 

[ y2(S(t-s)Q 1 / 2 e i ,e j fds 
J ° itj >i 

E E ( f s2 ^ - s ) ds ) ^ 1/2e - e tf ds 

i>l <>1 V» / 



< coEE^W 72 ^-) 2 

= C ||A-V(^)qV2||2 2(h)) 

where we have used ParsevaPs identity, (|2.13[) and (|2.16[) . By (|2.9[) we have that 
— l/p — n< —a, and thus using also (|2.4j) . 

\\A" 1/{2p) Q 1/2 \\ 2 C2{H) = TriA-ypQ) = Tr(A~ 1 / p ~ K A K Q) 

< Tr^-VP-^H^QIII^ < Tr(yl- Q )||^Q|| 6(ff) . 

Finally, the proof of the Holder regularity in time of u uses similar techniques and 
is omitted. □ 



Remark 2.8. Note that assumptions (|2.8[) — (|2.9[) arc stronger than the minimal as- 
sumption || A~ 1 ^ 2f ^Q 1 ^ 2 \\c2^H) < +°° needed for the existence of a mean squared 
continuous solution. One can replace (|2.8[) - (|2.9|) by 



\\A^' 2 Q^\\ c , iH) < +oo 

for some s > as a single main assumption on A and Q and obtain Holder regularity 
of order less than min(|, ^). 

3. Space discretization 



In this section we discretize ()2.1|) in space by a standard piecewise continuous 
finite element method. We refer to the monograph [19] for further details on fi- 
nite elements. We shall derive strong error estimates for the spatially semidiscrete 
problem for smooth initial data only imposing Assumption [T] on b. We will see later 
that for time discretization and also for the fully discrete scheme, we have to put 
further restrictions on b. Let {Th}o<h<i denote a family of triangulations of T>, with 
mesh size h > and consider finite element spaces {Vh}o</i<ij where VJ, C Hq(T>) 
consists of continuous piecewise linear functions vanishing at the boundary of T>. 
In order to derive the finite element formulation we look for a -valued process Uh 
such that 

J (du h (t), X )+ I b(t-s)(Vu h (t),V X )dsdt = (dW Q (t), X ), X&V h: t>0, 
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We introduce the "discrete Laplacian" 

(3.1) A h :V h ^V h , (A h ip,x) = m,W X ), 4>,xeV h , 
and the orthogonal projector 

P h :H^ V h , (P h f, x) = (/, X), X e V h . 

It is clear that the operator Ah is a positive definite bounded operator on Vh- Let 
us note also that using the definition (|3.ip of Ah , the following uniform inequality 
can be easily derived 

(3.2) \\A- 1/2 P h x\\ < WA-^xl x€H. 

Then, using the instability of Ph and some interpolation theory, we also have that 

(3.3) \\A- s P h x\\<\\A- s xl £e[0,i], xeH. 

Similarly to — A, the operator — Ah generates an analytic contraction semigroup on 
Vh and satisfies the uniform resolvent estimate 

\\z{z + Ah)' 1 P h \\ = \\zR(z,A h )P h \\ < 

for z G £</, = {z e C : |arg(z)| < <f> < 7r}. Since A h R(z,A h ) = I - zR(z,A h ), it 
follows that 

(3.4) \\AhR(z,Ah)Ph\\ B (H) < M<p + 1, ze S . 

Then we can rewrite the spatially semidiscrete problem in the same form as the 
original one as 



(3.5) 



du h + (J^ b(t - s)A h u h (s) ds^j dt = P h dW Q (t), t > 0, 
u h (0) = PhU . 

Similarly to the original problem the weak solution is given by 

u h {t) = Sh (t)P h u Q + f S h (t - s)P h dW Q (s), 
Jo 

where the resolvent family {Sh{t)}t>o can be written explicitly as 

oo 

Sh(t)P h u = y^^s h ,k{t){uQ,e ht k)e h ,k- 

k=l 

Here (\h.k,eh.k) are the eigenpairs of Ah and Sh,k(t) are the solution of the ODEs 
Sh,k(t)+^h,k / b(t - s)s h ,k(s)ds = 0, s/, fc (0) = l. 



We have the following stability result. 

Lemma 3.1. If b satisfies Assumption^ then for some C > 0, 



t 

2 J„ / /nrii ||2 



||S'(s)a;|| 2 ds < C\\x\\i L , t > 0, 

P 

\\S h {s)P h xf ds < C||i||ii, t>0, h>0. 
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Proof. We have, by ([2J3]l and (j2l6| . that 

/ \\S(s)x\\ 2 d S = J2 / s 2 (s)ds(x,e k ) 2 
Jo k=1 Jo 



< ^2\\ s k\\L^(R + )\\sk\\LHR + )(x,e k ) 2 < C ^2X k 1/p {x : e k ) 2 = C \\x\\ 2 _ L . 

k=l k=l 

As the constants in (|2.13|) and (|2.16|) do not depend on X k , we similarly obtain 



[ \\S h (s)P h x\\ 2 ds < C \\A- 1/2p P h x\\ 
Jo 



Finally, since -1/2 < -1/2,0 < -1/4, using ([53)1 with S = 1/(2/9), completes the 
proof. □ 

The error analysis is based on the Ritz projection 

R h : Hl(V) V h , (VR h v,V X ) = (Vu, V X ), w e x € 

In particular, we assume that Rh satisfies the error bound 

(3.6) \\R h v - v|| < C/i 7 ||w|| 7 , v g L>(A 7 / 2 ), 1 < 7 < 2. 

This puts some restriction on the domain T> but it is satisfied for convex polygonal 
domains, for instance. 

Next we prove an L 2 ((0, 00), H) error estimate for the space semidiscretization 
of the deterministic problem. It is an extension of the result in [2] where the special 
kernel b(t) = r^y e_ *' : ' 3_1 was considered. 

Proposition 3.2. If b satisfies Assumption^ and (|3.6[) holds, then 
\\S(t)x - S h (t)P h x\\ 2 dt < Ch 2s \\x\\ 2 _ k , < s < 2. 
Proof. It follows from Lemma 13. II that 

(3.7) / \\(S(s)-S h (s)P h )x\\ 2 ds<2 [ \\S(s)x\\ 2 + \\S h (s)P h x\\ 2 ds<C\\x\\\. 



Jo Jo 
To prove an error estimate of optimal order we set 

e(t) := S(t)x - S h {t)P h x := v{t) - v h {t) 



v(t) - P h v{t) + P h v{t) - v h {t) := p{t) + 9{t). 



For p, using the best approximation property of Ph, we obtain by Lemma 13. II and 
(13~B that 



(3.8) / \\ P (t)\\ 2 dt< / \\(R h -I)v(t)\\ 2 dt<Ch*\\x\\i_i. 

Jo Jo " 

In a standard way one derives an equation for 8 which reads 

6(t)+ f b(t-s)A h 9(s)ds = A h P h f b(t-s)(R h -I)v{s)ds, t > 0, 

Jo Jo 
(9(0) = 0. 

Taking Laplace transforms of both sides yields 

z6(z) + b(z)A h 6(z) = A h P h {R h - I)v{z)b{z). 



12 



MIHALY KOVACS AND JACQUES PRINTEMS 



Therefore, 

(3.9) 6{z) = A h R(^-,A h )P h (R h - I)v(z). 

b(z) 

It can be shown that b extends continuously to iM. \ {0}, see, for example, [14] . 
Therefore, using (|2"1)|) . it follows that ^ £ S , k £ R\ {0}, with <f> < it. Thus, 

\\A h R{^ y A h )P h \\ B(H) < (M + l)by (El. Therefore, setting 2 = **, k £ M\{0}, 

in (|3.9[) and using the isometry property of the Fourier transform we obtain, by 
Lemma 1331 and (|3.6|) . that 



(3.10) / ||#(i)|| 2 di < (A^ + l) / \\(R h - I)v(t)\\ 2 dt < Ch^xUi. 

Jo Jo p 

Interpolation using (|3 . T[) . (|3.8[) . and (|3.10[) yields 

/>oo 

|je(0|| 2 dt < 2 / (\\ P (t)\\ 2 + \\0{t)\\ 2 ) dt < C/i 2s ||ar||;_i, < s < 2. 



□ 

Next, using the error analysis from |13| . we have the following pointwise smooth 
data estimate for the spatially semidiscrete scheme. 

Proposition 3.3. Ifb satisfies Assumption^ and (|3.6I) holds, then for every e > 
there is C = C(T, e) such that 

\\S(t)x-S h (t)P h x\\ < C^ s |M| s(1+e) , < s < 2, t £ [0,T]. 

Proof. As already observed, Assumption[T]implies that b is a positive definite kernel. 
Therefore by, Q21 Theorem 2.1], it follows that 

/ rt 

,2 



\\S(t)x - Sh(t)Phx\\ < Ch \ \\x\\2 + J ||5(s)x||ada 
Proposition 12. 51 implies that 

(3.11) f \\S(s)x\\ 2 ds = [ \\A-*S(s)A 1+ *x\\ ds < C(T, e)\\x\\ 2+2t . 

Jo Jo 

Finally, since \\S(t) — Sh(t)Ph\\j3(H) < 2, interpolation finishes the proof. □ 

Theorem 3.4. Let A and Q satisfy (2.8\) - ll2.9\) and let b be satisfy Assumption^ 
//E||uo|j 2 ( - 1+£ - | < oo and (|3.6j) holds, then there is C = C(T,e,v) such that 

(E\\u(t)-u h (t)\\ 2 ) 1/2 <Ch v , u<--a + K, t£ [0,21 

9 

Proof. By the variation of constants formula, 

u(t) - u h (t) = S(t)x - S h (t)x + f (S(t-s)-S h (t-s)P h )dW Q (s). 



Thus, 

E||u(t) - u h (t)\\ 2 < 2E\\S(t)x - S h (t)xf 



2E 



{S(t-s)-S h (t-s)P h )dW Q (s) 



2 



ei + e 2 . 
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ei <Ch 2 »E\\uo\\l (1+er 



To bound e-i we use Ito's Isomctry and Proposition 13 . 21 to obtain 



(3.12) 



<?2 



(S(t-s)-S h (t-s)P h )dW Q (s) 



\\(S(t-s)-S h (t- S )P h )QV 2 \\l 2{H) ds 

= E / \\(S(s)-S h (s)P h )Q 1 ^e k \\ 2 ds 
k=i Jo 

OO 

< Ch 2 » \\A^' 2 Q x ^e k f = C/> 2 1^-^ 2 Q 1/2 !lLw 

k=l 

= Ch 2l, Tv(A u - L pQ) < Ch 2 "TY(A"- L P - K )\\A K Q\\. 



□ 



Remark 3.5. In particular, if Q = I, then d = 1, K = and a > \ whence v < - — ^. 
Also note, that it is clear from the proof that instead of (|2.8[) - (|2.9[) we could assume 
that ||^4 * Q^ll^m) < 00 ancl get a convergence rate of order v. Then, for 
trace class noise; that is, when Tr(Q) < oo we can take v = -. 



We end this section by showing that the above error estimate is optimal in the 
sense that it corresponds to the spatial regularity of the solution. 

Theorem 3.6. Let A and Q satisfy ifO)) - ifO)) and let v = - — a + k, or, let 

\\A^"~p^ 2 Q 1 / 2 \\c 2 ^h) < oo for some v > 0. If b satisfies Assumption^ and 
E\\u \\l < oo, then E\\u(t)\\l < C for some C > for all t > 0. 

Proof. It follows by Ito's Isometry and the fact that ||S(i)|| < 1 that 



u(t)\\l<2E\\u \\l + 2 / \\A^ 2 S(s)Q 1 / 2 \\ 2 :2iH) ds. 
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Let (efc,Afc) be the eigenpairs of A. Then, by monotone convergence, the self- 
adjointness of A and S, and Proposition 12. A\ it follows that 

/ j|^/ 2 s( s )Q^||2 ds = W n^wg^H^ 

Jo k=1 Jo 

OO 4 OO f 

J2 (A»/ 2 S(s)Q 1 / 2 e k ,e j ) 2 ds=Y^ [ (Q 1/2 e k , S^A^erf ds 
,k=i Ja j,k=i Jo 

oo ( t 

J2(Q 1/2 e k ,X; /2 e^ / s 2 (s)ds 
i_i Jo 



j 

oo 



OO 



< ]T (Q 1/2 e fe ,A; /2 ej ) 2 || Sj ||L^( R+ )||^||L 1(R+ ) 



"/2 \2i 



= Coll^-^^Q^ii^ < C Tr(^-^")||^Q||| (H) . 



□ 



4. Time discretization 



Time discretization is achieved via a classical implicit Euler scheme and, con- 
cerning the convolution in time, via a quadrature rule based on (|1.7[) . Let At > 
and we set t„ = n At for any integer n > 0. We seek for an approximation u„ of 
u(t n ) defined by the recurrence 

(4.1) Un-Un-X + At (^u n - k Au k \ = W Q {t n )- W(i n _i), n > 1, 

with initial condition uq = u(Q). We recall that the coefficients {c>Jk}k>o of the 
quadrature are chosen such that 

+oo /I _ \ 

(4.2) J2^z k = b(—^\, \z\<l. 

k=0 ^ ' 

Let us note that thanks to p~4j estimate (3.6)], we have the lower bound for ujq: 

(4.3) uj = b(l/At) > cAt p -\ At < 1, 
where p € (1,2) is defined in 



In the sequel we derive a discrete mild formulation (variation of constants for- 
mula) for (|4.ip . This formulation can not be made easily explicit as a function of 
the operators A, Q and the kernel b, because of the memory effect in the drift. First 
consider the deterministic algorithm 

(4.4) v n - v n _x + At cj n _ fc Avi^j = 0, n>l; v = x. 

Taking the z-transform, using the notation 

oo oo 

V{z) = VkZ h and ui(z) = ^kZ k , 

k=0 k=0 
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we get 

V(z) - x - zV(z) + Atu(z)A(V(z) - x) = 0. 

Thus, 

V(z) = (I + Atui(z)A)((l - z)I + Attu(z)A)~ 1 x := B(z)x, 

where 

oo 

B(z)x = ^B k xz k . 

k=0 

This means that Vk = Bkx, k = 0, 1, ... Note that Bq = B(0) = I. For the stochastic 

equation it will be useful to rewrite B(z)x as 

(4.5) 

B(z)x =((1 - z)I + u>(z)AtA)~ 1 (I + Lo(z)AtA)x 

= ((1 - z)I + Lj(z)AtA)- 1 x + u(z)AtA((l - z)I + Cj(z)AtA)- 1 x 
= ((1 - z)I + uj(z)AtAy 1 x - (1 - z)((l - z)I + u{z)AtAy 1 x + x 
= (z((l - z)I + uj(z)AtA)- 1 + I)x. 

Now, we consider the stochastic case (|4.ip which reads, after taking the z-transform, 
rearranging, and using the notation w„ = W^(t n ) — W^{t n -\) for n > 1, Wq = 0, 
and 

oo oo 

w{z) = 2_, w kZ n and U(z) = V] UkZ k , 
k=a k=a 

as 

U(z) = B(z)x + ((1 - z)I + uj(z)AtA)- 1 w(z) 

= B(z)x + B< " Z ^ I w(z) = B(z)x + B{z)^- - -w{z), 

where we also used (|4.5[) to rewrite the stochastic term in the previous calculation. 
This yields the discrete variation of constants formula, taking into account that 
wq = and that Bq = I, 

n n—1 
(4.6) U n = B n X + ^ Bn-kWk+l - W n+ i = B n X + ^ B n-kWk+l- 

k=0 k=0 

The importance of this formula lies in the fact that it connects the deterministic 
case to the stochastic case with the deterministic time-discrete solution operators 
B n explicitly appearing in the formula. 

4.1. Deterministic estimates: stability and smoothing. The next theorem is 
interesting in its own right. It shows that Lubich's convolution quadrature based on 
the Backward Eulcr scheme have a remarkable qualitative property: it preserves the 
L p -norm of the orbits of the solution. The result can be viewed as a generalization of 
the ones in [6]; in particular, it removes the additional technical frequency condition 
in [6l Theorem 2] . The proof uses a representation similar to that in [I] . We also 
note that the statement holds in Banach spaces as well since the proof does not use 
Hilbert space techniques. 
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Theorem 4.1. If the resolvent family {S(t)}t>o of ()2.10j) satisfies 

S(-)x g L p {{0,oo):H) 
for some 1 < p < oo and x € H , then 



and 



At V \\B k x\\ p < / HS^xfdt, l<p<oo, 

fc=l Jo 



sup||.B fe x|| < ||5(-)x|| iro(R ). 

k>l 



Proof. The Laplace Transform of {S(t)}t>o is given by 

S(z)x = (zl + b(z)A)- 1 x. 
Using (|4.2[) and (|4.5|) we see that the z-transform Bx of {B n x} n is given by 

1 - 1 — z f°° 
B(z) = z—S(——)x +x = x + z / S(Ats)e- s e zs ds 
At At Jq 

Therefore, we conclude that Bq = I and that 

(4.7) B k x = f S(Ats)x e ,, — ,| ds for fe > 1. 

Jo 0-1)! 



Let 



(fc-1)! 

Then > 0, ||AHl 1 (r + ) — 1- Therefore, if p = oo, we immediately obtain from 

mm that 

sup||.B fc x|| < \\S(-)x\\ L ^( R+) . 

k>l 

If 1 < p < oo, then we use Jensen's inequality in (|4.7I) . and have 

n n „oo 

At^||£ fc x|| p <]TAt / ll^(A^ll p /fe(s)d S 
fc=i fe=i 70 

= ||5(t>r^/ fe (— )dt< B up£/*(t) / ii5(t>p<tt. 

Jo fc=1 At t>0 ^ J 

Finally, noticing that X)^Li /« = 1 completes the proof. □ 

Theorem 14. II has the following important corollary on the smoothing and stability 
of the time discretization scheme in case b satisfies Assumption 1. 

Corollary 4.2. If b satisfies Assumption^ then, for all x € H, 

n 

sup||-B fc x|| < ||x|| and At V ||-B fc x|j 2 < C||x||ii, n > 1. 
fe ^ ft 

Proof. The statement follows from Theorem 14. II together with Lemma |3. II and the 
fact that || S(t)\\ < 1 for t > 0. □ 



Finally we will need a Holder type estimate on the resolvent family {S'(t)} 



t>o- 
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Lemma 4.3. If b satisfies Assumption^ then there is C — C(T, 7) > such that 
/"/■** \ 1/2 

(J2 \\(S(t n -s)-S(t n -t k . 1 ))xfds\ < CA^||x|| s _i, nM = T, 
for all 7 < 4^ where < s < - . 

Proof. It follows from (|2.17p . with s = and Lemma lXTl that there is a constant 
C = C(e, T) such that, for < e < ±, 

(n t \ 

X) / ||(^„- S )-5(i„-t fc -i))x|| 2 d S ) <C||a;|| e _i,nAt = t„=r. 
fe=1 ->tfe-l j 

Next, it follows from Proposition 12 .41 that 

/ ||(5(^- S )-5(t„-t fc _i))x|| 2 ds 
fc=1 "'tfc-i 

00 n tk 

^2(x, e-i) 2 ^ / " s) - s 4 (t„ - i fc _i) 2 ds 

^ 1 1 i Jtu_-\ 



ri 



<2^(x, ei) 2 ^ / |sj(t„ - s) - Sj(t„ - *fe_i ) I ds 
i=i fe=i^*fc-i 

00 n (-tfc rt n -t k -l 

< 2j2( x ^ifJ2 / / |si(t)|dtda 

-■1 r 1 J tn— S 



i=l k=l* ,I »'-i 



fe=l 



tfc ,.t„-t 



<2^(x,e 4 ) 2 ^ / / [s 4 (t)|d*da 



A' - 1 



tfc— 1 Jtn—tk 



< 2At||x|| 2 sup||s I || L i (R+) < CAt\\x\\ 2 . 

Finally, interpolation gives the desired result. □ 

4.2. Deterministic estimates: convergence rates. In order to give an error 
estimate of optimal order with no initial regularity for the time discretization of 
the deterministic problem we have to impose another assumption on b. This kind of 
assumption; that is, the existence of an analytic extension of & to a sector beyond 
the left halfplane, is fairly standard in the existing deterministic literature, see, 
for example, [5J [HI [TTJ [T^], but it clearly represents a major restriction compared 
to Assumption [TJ We note that this additional assumption is not needed neither 
for the spatial error estimates with smooth initial data, and hence for the spacc- 
scmidiscrctization of the stochastic equation, nor the for the stability results for 
the time discretization in the previous subsection. 

Assumption 2. The Laplace transform bofb can be extended to an analytic func- 
tion in a sector S e with 6 > tt/2 and |S (fc) (z)| < C\z\ 1 ^ p ' k , k = 0,1, z £ Eg. 

Note that Assumption [5] implies that 

(4.8) w Q = 6(1/ At) < CAt p -\ At < 1. 

An important example of a family of kernels satisfying both Assumptions [T] and [2] 
is given by 6(t) = Ct^e^, < /? < 1 and 77 > 0. 



18 MIHALY KOVACS AND JACQUES PRINTEMS 

Assumptions Q] and [2] allows us to use the following deterministic nonsmooth 
data estimate [TSJ Theorem 3.2]. 

Proposition 4.4. If Assumptions^ and& hold, then there exists C = C(p) > 
such that 

(4.9) \\S(t n )x-B n x\\<^At\\xl n>l. 

Corollary 4.5. If Assumptions^ and&hold, then there exists C = C(T, 7, p) such 
that 

/ \ 1/2 

At]T \\S(t k )x -B k x\\ 2 \ <CAF\\x\\ B _i, nAt = T, 

V k=0 ) 

for all 7 < inhere < s < - . 

Proof. It follows from (|2.17j) . with s = - — e, and Corollary 14.21 that there is a 
constant C = C(e,T) such that, for < e < -, 

/ \ 1/2 

f AtJ2j\S(tk)x-B k xf) < C\\x\\ e _L, nAt = T, e > 0, 

where we also used the fact that Bq = 5 (to) = I- Furthermore, since \\S(t k ) — B k \\ < 
2 by Corollary K2\ it follows from (gjj]) that 

||5(t fc )a? - B k x\\ <CAtt-H e -*\\x\\, k>\, 
and thus, for some C = C(e, T, p), 

/ „ x 1/2 

AtJ^ ||S(t fc )i - B fc x|| 2 < CAti- e \\x\\. 

V fc=o / 
Interpolation finishes the proof. □ 

4.3. Error estimate for the stochastic equation. We can now state and proof 
the main result of this section. 

Theorem 4.6. Let A and Q satisfy (|2.8[) - (|2.9p and let b satisfy Assumptions [7] 
and\M Suppose further that E||uo|| 2 < 00. For T > 0, let {u(t)}te[o,T] be the 
unique weak solution of (|2 . 1 [) and let u n be the solution of the scheme (|4.1j) with 
T = nAt. Then for any 7 < (1 — p(a — k))/2, there is C = C{p, E|jito|| 2 ) > and 
K = K{T, a, 7, k, p) > such that 

(4.10) (E|KT) - u n \\ 2 ) 1/2 < CT~ l At + KAt 1 1 t„ = nAt = T. 
Proof. If e n = u(T) - u n = u(t n ) - u n , then ([2331) and $M yields 



e„ = (S(t n ) - B n )u + 2J 



tfc 

(S(t n -s)-B n _ k+1 )dW°(s) 

t k -i 



k=l L 

Taking the expectation of the square of the .ff-norm of e n leads to, by independence 
and Ito's isometry: 

(4.11) E||e„|| 2 < 2(a + 6), 
where a denotes the deterministic part of the error: 

(4.12) a = E\\(S(t n )-B n )u Q \\ 2 , 
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1!) 



and b the stochastic part 

+00 n 



b = J2J2 \\(S(t n -s)-B n _ k+1 )Q^ 2 e i \\ 2 ds. 

i=l fc=l 



Thanks to (|4.9p . a can be bounded as 

(4.13) a < ■^At 2 E||u || 2 , n > 1. 

We use Corollary 14 . 5 1 and Lemma T4. 31 to bound b as 

b<^J212 \\(S(t n -s)-S(t n -t k ^ 1 ))Q 1 / 2 e i \\ 2 ds 
i=i fe=i Jtk-i 

00 n „t fc 

+ \\(S(t n -t k .. 1 )-B n - k+1 )Q 1 ' 2 e i \\ 2 ds 

OO 

< CAt* Y, \\Q 1/2 ei\\ s -i = CAt^\\A^-^ 2 Q^Yc 2W 

i=i 

< CA^Tr(^ s -^- K )||A K Q||2 (H) . 

Finally, we set —a = s — - — K and conclude that 7 < ^ = (1 — p(a — n))/2. □ 



Remark 4.7. In particular, if Q = I then d = 1, ft = and a > \ whence 
7 < 1/2 — |. Also note, that it is clear from the proof that instead of (|2.8j) ~ ()2.9j) 

we could assume that ^A^ ~~^ 2 Q 1 ^ 2 \\c->(H) < °o an d obtain 7 < Then, for 
trace class noise; that is, when Tr(Q) < 00 we can take s = - and hence 7 < 1/2. 
Remarkably, this is the same rate as for the heat equation |21j independently of the 
value of p. 

5. The fully discrete scheme 

In this section we derive strong error estimates for a fully discrete scheme for 
(|2.1[) . Both Assumptions 1 and 2 on b are needed but in return we get optimal 
error bounds with no initial regularity. As the fully discrete scheme, similarly to 
the time semidiscretization (|4.ip . we consider the recurrence 

(5.1) Vft-^-i.h + A^ ^wn-tipuj = Ph(W Q (t n )-W c ^(t n - 1 )), n>l, 

with uo,h = PhUQ. Again, the solution is given by the discrete variation of constants 
formula 

ra-l 

(5.2) u n , h = B n:h P h u + Bn-^hPhAW^, 

fc=0 

where AW k +i = W(t k +i) ~ W{t k ) and {B kt h}k>o is a family of linear bounded 
operators with Bq ^ = I. 

Theorem 5.1. Let A and Q satisfy H2.8\) ~ nT3\) and let b be satisfy Assumptions^ 
and\^ Suppose further that E||uo|| 2 < 00. ForT > 0, let {u(t)}te[o,T] be the unique 
weak solution of (|2.1[) and let u n _h be the solution of the scheme (|5.1[) with T = nAt. 
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// (3J| hol ds, then there is C = C(p,E||u || 2 ) > and K = K(T,a,-r,K,p) > 
such that 

(5.3) (E||u(T) - u n>h \\ 2 f' 2 < CiAtT- 1 + h 2 T- p ) + K(Af< + h u ), nAt = T, 
where 7 < (1 — p(a — n))/2 and v < - — a + n. 
Proof. We decompose the error as 

u(T) - u n< h = S(T)u - B n . h P h u 



f S{T-s)dW Q (s)- f S h {T-s)P h dW Q {s) 

Jo Jo 

f S h (T - s)P h dW Q (s) - V B n - kih P h ,AW fe ° +1 

JO ,. r, 



k=0 

:= ei + e 2 + e 3 . 

First we bound e± which is the deterministic error. Under Assumptions [JJ and [5] we 
have that 

(E|| ei || 2 ) 1/2 < C{AtT- 1 + tfT-P-^iEWuof) 1 ' 2 
by [T^J Theorems 2.1 and 3.2]. Next, e 2 has already been bounded in (|3.12p as 

(5.4) E||e 2 || 2 < Ch 2 "\\A (u -^ /2 Q 1/2 \\ 2 C2(H) < Ch 2l/ Tr(A l/ ~'p~ K )\\A K Q\\. 

Finally, the proof of Theorem 14.61 shows that, 

(5-5) E|| e3 || 2 < CAt^\\A { ^ ), \p h QP h f'Yc 2{H y 

Set —r = (s — ~)/2 and note that since < s < - we have that < r < 1/2. Then, 

\\A- h r {P h QP h ) 1/2 \\l 2{H) = T,{P h A- h r P h QP h A;SP h ) = WA^P^Ymh) 
<\\AZ r P h A r \\ 2 B{H A\A- r Q^ 2 \\ 2 C2{H} . 
Thanks to (O with 6 = r € [0, 1/2), it follows that \\A~ r P h A r \\ B{H) < 1. Hence, 

E|| e3 || 2 < CAt^WA^-^Q^fc^H) < CAt 2 ^{A s -T- K )\\A*Q\\ 2 B(HV 
and the proof is complete. □ 

Remark 5.2. We would like to highlight two important special cases. Firstly, if 
Q = I then d = 1, K = and a > \. Hence v < - - \ and 7 < 1/2 - |. As 

before, we could assume, that \\A {v -^ l2 Q l / 2 \\ c ^ H) < 00 instead of $nj$ and (1231) 
and get a convergence rate of order v is space and 7 < ^ in time. In particular, if 
Tr(Q) < 00, then we may set v = -. Thus, the time order is almost 1/2, the same 
as for the heat equation with trace class noise, but the space order is less than 1, 
which is the space order for the heat equation, see j21) . 

Remark 5.3. The pure time-discretization as well as the fully discrete scheme can 
be studied for smooth initial data under Assumptions [1] and [2] on b. Using (TBI 
Theorem 3.1] and |12[ Lemma 3.2] one arrives at the deterministic estimate 

(5.6) \\S(T)u Q -B nth P h u \\ 

<C(T)(h 2 + k)(\\u \\ 2 + [ T \\S( S )u a \\ 2 ds+ [ T \\S(s)u \\ds ) . 
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If ito <E 'D(A), then u{t) = S(t)uo is a strong solution of (|2.10j) . see JXT1 Proposition 
1.2]; that is, u(t) = S(t)u satisfies (f2~T0| with / = for all t > 0. Then 

S(t)u + [ b(t- s)AS(s)u ds + b(t)Au = 0, t > 0, 
Jo 

and thus 

jf ||£(«)u || ds < C(T) \\S(s)u \\ 2 ds+\\u \\ 2 Y 

Therefore, using stability, interpolation and (|3.1ip . it follows that 

\\S{T)u - B n , h P h u \\ < C(T,e)(h s + k s ^)\\u \\ sil+e)l 0<s<2. 

The latter estimate can be used to replace the first term in the bound (15.31) in case 
E||uo|!s(i +e ) < 00 ' < s < 2. The estimates for the pure time-discretization arc 
analogous using [12l Theorem 3.1] which states (|5.6|) with h = and B n> o = B n . 
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